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ABSTRACT 

Excitation of a Global Alfven wave (GAW) is proposed as a viable mechanism 
to explain plasma heating in the magnetosphere of young stars. The wave and 
basic plasma parameters are compatible with the requirement that the dissipation 
length of GAWs be comparable to the distance between the shocked region at 
the star's surface and the truncation region in the accretion disk. A two-fluid 
magnetohydrodynamic plasma model is used in the analysis. A current carrying 
filament along magnetic field lines acts as a waveguide for the GAW. The current 
in the filament is driven by plasma waves along the magnetic field lines and/or 
by plasma crossing magnetic field lines in the truncated region of the disk of the 
accreting plasma. The conversion of a small fraction of the kinetic energy into 
GAW energy is sufficient to heat the plasma filament to observed temperatures. 

Subject headings: accretion, accretion disks — magnetized plasmas — electron- 
ion collisions — Alfven Waves — MHD eigenmodes — driven current — young 
stars 



1. Introduction 

In the present popular model for classical T Tauri stars (CTTSs), the so called class II 
objects (Lada 1987), a central young star is surrounded by a geometrically thin accretion 
disk. The disk is disrupted at a given radius by the magnetic field of the star. For smaller 
radii the accretion flow follows the magnetic field lines of the star until it impacts onto 
the stellar surface Ghosh & Lamb (1979); Konigl (1991). This magnetospheric accretion 
model explains observational signatures seen in some CTTSs as, for example, the excess 
of optical and ultraviolet continuum flux (veiling), as well as redshift absorption features 
in emission line profiles (MuzeroUe et al. 1998; Hartmann, Hewett & Calvet 1994). The 
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plasma of the accretion disk flows along the dipole magnetic field lines towards the magnetic 
poles of the star. It is inhomogeneous in the plane perpendicular to the magnetic field 
lines, but although somewhat more homogeneous along the field lines, cavities and plasma 
condensations (plasmoids), aligned by the magnetic field, may be present. These plasmoids 
are distributed randomly in the plane perpendicular to the plasma flow. Waves excited by 
the shock of the plasma flux with the star's surface can propagate upstream and interact 
with the plasmoids, creating a current, which in turn, will create current - carrying filaments. 
An upstream or downstream electric current can be driven by the waves or by the v x B- 
MHD effect in the truncated region of the disk. Such filaments have been observed in the 
laboratory (see, e.g. Fadeev et al. (1965); Rosenberg et al. (2001)). The occurrence of 
filaments has also been predicted in simulations of rotating magnetized disks (Machida et al. 
2000). 

Hartmann, Hewett & Calvet (1994) and, more recently, Muzerrole et al. (1998) carried 
out calculations on magnctosphcre accretion models, solving the radiative transfer equations 
in the Sobolev approximation. Theoretical Balmer line profiles obtained were found to be in 
reasonable agreement with observations. The major uncertainty in radiative magnetosphere 
model calculations is the temperature profile along the tube. Martin (1996) calculated the 
temperature profile from the energy balance of the gas, including heating by adiabatic com- 
pression of magnetic field lines, Balmer photoionization, and ambipolar diffusion. However, 
the values of the temperature that he obtained were too low to explain the observed spectra. 
Thus, additional heating mechanisms inside the tube are needed. 

The line-averaged electron temperature in the magnetosphere taken from optical mea- 
surements, is about 0.8 eV (see, e.g. Hartmann, Hewett & Calvet 1994; Martin 1996; Muzer- 
role et al. 1998). Balancing the heating rate with the black body emission at a temperature 
^ 0.8 eV, we find that about 5% of the fiux energy is sufficient to heat the downstream fiux 
to a few eV. The main consumption of the flux energy is spent on plasma ionization. Us- 
ing laboratory experimental data for the ionization - recombination processes (see Delcroix 
1965), we can estimate the energy necessary for the plasma ionization rate to be about 50% 
in the fllamcnt for electron temperatures of about 1.6 cV. We use here a theoretical tem- 
perature proflle, varying from ~ 7500 K to ~ 8300 K, to reproduce the optical line features 
observed. 

Taking into account the evidence for turbulence in the observations, as well as the 
necessity for an additional heating mechanism, Vasconcelos, Jatenco-Pereira & Opher (2002) 
suggested that Alfven waves may be important in the heating of magnetic flux tubes of 
CTTSs. They studied the possibihty that the waves generated at the star's surface due to 
the shock are produced by the accreting matter, as suggested by Scheurwater & Kuijpers 
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(1988). 

Various damping mechanisms for Alfven waves have been suggested in the hterature, 
such as Alfven resonant heating of solar loops (HoUweg 1990; Elfimov et al. 1996), wave 
damping by phase mixing (also in the solar context), and cyclotron heating, occurring as an 
Alfven wave travels down a magnetic field gradient until its frequency matches that of the 
cyclotron resonance of helium or some other plasma specie due to the decreasing ion-cyclotron 
resonance (magnetic beach). In addition to the more conventional collisional and viscous- 
resistive Alfven wave damping (Osterbrock 1961), Vasconcelos et al. (2002) concentrate on 
nonlinear and turbulent damping (Vasconcelos, Jatenco-Pererira & Opher 2000). Nonlinear 
and turbulent damping have also been studied by Jatenco-Pereira & Opher (1989a); Jatenco- 
Pereira, Opher & Yamamoto (1994) in the solar wind, Jatenco-Pereira & Opher (1989b) in 
protostellar winds, Jatenco-Pereira & Opher (1989c) in late-type giant stars, dos Santos, 
Jatenco-Pereira & Opher (1993a; 1993b) in Wolf-Rayet stars, Gongalves, Jatenco-Pereira & 
Opher (1993a; 1996) in quasar clouds and Gongalves et al. (1993b) in extragalactic jets . 

The (GAW) has not been previously analyzed as a possible mechanism for the energy 
transport from the shock and truncation regions to the magnetosphere of T Tauri stars. 
They were predicted in numerical calculations (Ross et al. 1982) using a cylindrical model 
for magnetized plasmas with an axial current and were observed later in tokamak experi- 
ments. GAWs are excited in plasmas in a helical magnetic field, with the poloidal component 
produced by an axial current, -Be.e/ ~ 27iaBz/L, where is the axial magnetic field com- 
ponent and a and L are the plasma radius and length, respectively. The so-called tokamak 
safety factor, q = rBz/RBg, is considered to be unity. In cold plasmas with temperatures 
less than ~ 0.1 eV, GAWs dissipate via electron- ion coUisions. However, in hot plasmas, 
when the Alfven velocity, ca — -B/(47rmine)^/^, is of the order of the electron thermal veloc- 
ity, VTe — {Te/meY^'^, GAWs dissipatc via electron Landau damping. The electron Landau 
damping mechanism is very effective when waves interact with a group of electrons that 
have velocities approximately equal to the wave phase velocity, uj/k\\ = ca, where u is the 
frequency and fcy is the wavenumber parallel to the magnetic field lines. Electromagnetic 
field components related to GAWs are proportional to exp[i{m9 + kzZ — ujt)], with poloidal m 
and axial wave numbers. The requirement for these waves to propagate is that m = — 1 
and < 0, in the case of current flowing along the magnetic field. A GAW has a very 
large quality factor, ratio of frequency to dissipation decrement, Q = cu/j 1 and a 
very long dissipation length in comparison to their wave length. La ^ Thus it is 

a very good candidate to explain plasma heating far from a source, such as in solar loops 
(HoUweg 1990; Elfimov et al. 1996). In addition a GAW is able to drive the current, which 
can be calculated by balancing the wave momentum transfer force with the electron-ion 
friction force (Elfimov et al. 1996). 
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2. Magnetosphere Accretion Model 



We study here a two-fluid magnetohydrodynamic model of a young star, accreting plasma 
(see Fig. 1). The standard parameters of CTTSs, similar to those used by Vasconcelos 
et al. (2002), are assumed: star radius = 10^^ cm, truncated radius Rfr = 2 x 10^^ 
cm, truncated depth A^^ ~ 3 x 10^° cm, magnetic field at the star's surface Bgt = 1 kG 
and Bfr ~ 0.1 kG in the truncated region, where the electron density is fitr — 3.0 x 10^^ 
cm~^. The length L^nf of the magnetic field hues from the shock at the star's surface to the 
truncated region is ~ 1.2R^. The velocity of the shock is Vsh ~ 2.5 x 10^ cm at the star's 
surface and Vtr ~ 3.5 x 10'' cm s~^ in the truncated region. The ionization rate parameter 

< av >ion for Tf. = 1.5 eV is ~ 3 x 10^^^ cm^ s~^ (Delcroix 1965). The accretion energy 
density is deflned by the ratio of the energy flux density to the path length L^f- Balancing 
ninety percent of the accretion power density minovfj./Lmf with the ionization power density 

< av >ion ngUoEion, we obtain a plasma with an average density along the magnetic field of 
n 3.3 X 10^^ cm-^ 

We assume that the magnetic fleld connecting the truncated disk and the star's surface 
consists of fllaments, along which the GAW propagates. The filament current can be driven 
by the wave ponderomotive force (see discussion below) and/or by a magnetodynamo in the 
truncated region of the disk. The current, I[a] — ^dB^ ef, and current density, joh — I /t^o? ~ 
16 statA cm^, in a filament of radius a can be produced by the electric field, j/a ^ 10~^^ 
statV cm"-*^, that is sufficient to create the required poloidal magnetic field B^^ef for the 
GAW excitation. This electric field can be generated by the induction of a perpendicular 
electric field, Ej_ = BtrVtr/c ~ 0.05 statV cm~^, in the truncated region of the disk due 
to the magnetodynamo v x B effect (see Fig. 1). A difference of potential on the scale of 
the filament radius is of the order of 2.5 x 10"^ statV. This potential is much higher than is 
necessary and can generate a filament current with a hot electron tail and a characteristic 
velocity higher than that of the Alfven velocity ca, as was found in the plasma focus device. 
We expect that the electrons in the tail will remain hot in the filament by their interaction 
with the GAW. 

For the parameters of the star assumed above, the phase velocity of a GAW, CA{Rtr), 
is ~ 3 X 10^ cm s~^ at the center of the filament. Taking into account the plasma resistance 
(see, e.g. Ginsbrug 1961), the dissipation length of the Alfven waves can be estimated to be 
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acA 



1.2 X 10^^ cm. 



(1) 
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1.6 X 10 
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= 3 X 10"^ cm-\ 



- 5 - 



This A; 1 1 is a maximum value, determined by the condition that the dissipation length be 
sufficient so as to explain the plasma heating along the magnetic field from the truncated 
region to the star's surface. 

Assuming a conversion of 5% of the total accreting plasma energy fiux, S — {l/2)nnmi{Vgf^+ 
vfj.), into the Alfven wave flux, Sa — Pa La, we obtain a wave dissipation rate of Pa — 0.0039 
erg cm~^ s^^. This dissipation rate produces a wave momentum transfer force PAk\\/uj (see 
Elfimov et al. 1996). Balancing the wave momentum transfer force with the electron-ion 
friction force UeimefieVe, we obtain the density of the driven current, 

< Jed >= — ■ (2) 

Using a reduced friction, fei,e/ ~ ^ei{vTe/ ca)^^^ , for the electron tail with effective velocity 
~ ca and u = k\\CA = 9 x 10^ s~^ {ca = 3. x 10^ cm s— 1), we obtain the required current 
density, < jcd >~ 16 statA cm~^, in the filament. In this case, the electron tail that 
transports the current along the filament can be supported by the GAW via the electron 
Landau damping, while the current in the filament can be driven by the GAW. 

Here, we also present the results of calculations for the GAW in the filament. In order 
to calculate the wave field and dissipation, we use the eikonal model to analyze the wave 
propagation along the filament, exp[i(m^^ + J k^dz — ujt)], where z is the coordinate along 
the filament axis and 6 is the poloidal angle. In this case, in order to simplify the GAW cal- 
culations using numerical codes (Amarante-Segundo et al. 1999; Galkin et al. 2002), instead 
of a dipole magnetic field, we assume a cylindrical model for the filament plasma with the 
current flowing along a magnetic fleld that has a helical magnetic fleld line conflguration. 
To calculate the wave characteristics along the fllament, we take the major radius Rq — 10^^ 
cm comparable to the star's radius and the minor radius a — 5 x 10^ cm. The calculations 
are carried out with the electron density profile of the form ng — no(l ~ jo?^ + ria, with 
no = 5 X 10^^ cm~^ and = 5 x 10^° cm~^. The temperatures are taken as Teo = 1.6 
eV and Tjo = 0.5 eV. The ion density is chosen so as to satisfy the requirement of charge 
neutrality, Znz + rii = rig, where rii and nz are the densities of hydrogen and the other 
ions, respectively. In accordance with the Spitzer conductivity oc Te^^, the current density 
profile has the form j — jo(l — r^/a^)^. The value of the safety factor, q — rBz/RBe, is 
assumed to be unity at the filament axis and to monotonically increase up to four at the 
plasma boundary. In Fig. 2, we show the characteristic profile of the magnetic field as well 
as the dissipated power of the GAW for the filament parameters discussed above. 

In Fig. 2, we observe the peaking of the dissipation at the filament center, which predicts 
filament center peaking of the electron temperature distribution, as well. The GAW heats 
the core of the filament to higher temperature, while reducing the plasma temperature at 



-6- 



the filament border due to heat diffusion. Tfie axial wavenumber of the GAW, kz — u/ca, 
varies slowly along the magnetic field lines since it depends on changes in the local magnetic 
field and density. 

3. Conclusions 

Our analysis of GAWs in the magnetospheres of CTTSs shows that: 

a) The energy of the shock (and/or accretion disk region) can be transported from the surface 
of the star (and/or accretion disk region) to the magnetosphere region by GAWs via current 

carrying filaments; 

b) The dissipation length of GAWs is on the order of the truncation radius; 

c) The conversion of only a few percent of the shock energy into GAWs is required in order 
to heat the filaments up to the observed temperature of the star's magnetosphere; 

d) The current in the filament may be induced by the electric field, created by the v x B- 
magnetodynamo effect in the truncated region of the disk. However, the major part of the 
current in the filaments is driven by GAWs. 

The authors thank the Brazihan agency FAPESP (Proc. No. 00/06770-2) for support. 
VJP and RO thank the Brazihan agency GNPq and the project PRONEX (41.96.0908.00) 
for partial support. 
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Fig. 1. — Sketch of a young accreting star of radius Ri, and truncation depth Atr- The 
electric field, designated by ®, is generated by the v x B-effect due to the interaction of the 
magnetic field lines from the star with the accreting plasma in the truncated region of the 
disk. 

Fig. 2. — Plot of Sj.^5)-components of the GAW field and dissipated power across the filament 
cross-section for a frequency of 1.4 kHz and the parameters B — IkG, go = 1-0, no — 5x 10^^ 
cm-3, Teo = 1.6 eV, and = 0.5 eV. 
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1. Introduction 

In the present popular model for classical T Tauri stars (CTTSs), the so called class II 
objects (Lada 1987), a central young star is surrounded by a geometrically thin accretion 
disk. The disk is disrupted at a given radius by the magnetic fleld of the star. For smaller 
radii the accretion flow follows the magnetic fleld lines of the star until it impacts onto 
the stellar surface Ghosh &; Lamb (1979); Konigl (1991). This magnetospheric accretion 
model explains observational signatures seen in some CTTSs as, for example, the excess 
of optical and ultraviolet continuum flux (veiling), as well as redshift absorption features 
in emission line profiles (MuzeroUe et al. 1998; Hartmann, Hewett & Calvet 1994). The 
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plasma of the accretion disk flows along the dipole magnetic field lines towards the magnetic 
poles of the star. It is inhomogeneous in the plane perpendicular to the magnetic field 
lines, but although somewhat more homogeneous along the field lines, cavities and plasma 
condensations (plasmoids), aligned by the magnetic field, may be present. These plasmoids 
are distributed randomly in the plane perpendicular to the plasma flow. Waves excited by 
the shock of the plasma flux with the star's surface can propagate upstream and interact 
with the plasmoids, creating a current, which in turn, will create current - carrying filaments. 
An upstream or downstream electric current can be driven by the waves or by the v x B- 
MHD effect in the truncated region of the disk. Such filaments have been observed in the 
laboratory (see, e.g. Fadeev et al. (1965); Rosenberg et al. (2001)). The occurrence of 
filaments has also been predicted in simulations of rotating magnetized disks (Machida et al. 
2000). 

Hartmann, Hewett & Calvet (1994) and, more recently, Muzerrole et al. (1998) carried 
out calculations on magnetosphere accretion models, solving the radiative transfer equations 
in the Sobolev approximation. Theoretical Balmer line profiles obtained were found to be in 
reasonable agreement with observations. The major uncertainty in radiative magnetosphere 
model calculations is the temperature profile along the tube. Martin (1996) calculated the 
temperature profile from the energy balance of the gas, including heating by adiabatic com- 
pression of magnetic field lines, Balmer photoionization, and ambipolar diffusion. However, 
the values of the temperature that he obtained were too low to explain the observed spectra. 
Thus, additional heating mechanisms inside the tube are needed. 

The line-averaged electron temperature in the magnetosphere taken from optical mea- 
surements, is about 0.8 eV (see, e.g. Hartmann, Hewett & Calvet 1994; Martin 1996; Muzer- 
role et al. 1998). Balancing the heating rate with the black body emission at a temperature 
^ 0.8 eV, we find that about 5% of the fiux energy is sufficient to heat the downstream fiux 
to a few eV. The main consumption of the flux energy is spent on plasma ionization. Us- 
ing laboratory experimental data for the ionization - recombination processes (see Delcroix 
1965), we can estimate the energy necessary for the plasma ionization rate to be about 50% 
in the fllament for electron temperatures of about 1.6 eV. We use here a theoretical tem- 
perature proflle, varying from ~ 7500 K to ~ 8300 K, to reproduce the optical line features 
observed. 

Taking into account the evidence for turbulence in the observations, as well as the 
necessity for an additional heating mechanism, Vasconcelos, Jatenco-Pereira & Opher (2002) 
suggested that Alfven waves may be important in the heating of magnetic flux tubes of 
CTTSs. They studied the possibihty that the waves generated at the star's surface due to 
the shock are produced by the accreting matter, as suggested by Scheurwater & Kuijpers 
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(1988). 

Various damping mechanisms for Alfven waves have been suggested in the hterature, 
such as Alfven resonant heating of solar loops (HoUweg 1990; Elfimov et al. 1996), wave 
damping by phase mixing (also in the solar context), and cyclotron heating, occurring as an 
Alfven wave travels down a magnetic field gradient until its frequency matches that of the 
cyclotron resonance of helium or some other plasma specie due to the decreasing ion-cyclotron 
resonance (magnetic beach). In addition to the more conventional collisional and viscous- 
resistive Alfven wave damping (Ostcrbrock 1961), Vasconcelos et al. (2002) concentrate on 
nonlinear and turbulent damping (Vasconcelos, Jatenco-Pcrcrira & Ophcr 2000). Nonlinear 
and turbulent damping have also been studied by Jatenco-Pcrcira & Opher (1989a); Jatenco- 
Pereira, Opher & Yamamoto (1994) in the solar wind, Jatenco-Pereira & Opher (1989b) in 
protostellar winds, Jatenco-Pereira & Opher (1989c) in late-type giant stars, dos Santos, 
Jatenco-Pereira & Opher (1993a; 1993b) in Wolf-Rayet stars, Gongalves, Jatenco-Pereira & 
Opher (1993a; 1996) in quasar clouds and Gongalves et al. (1993b) in extragalactic jets . 

The (GAW) has not been previously analyzed as a possible mechanism for the energy 
transport from the shock and truncation regions to the magnetosphere of T Tauri stars. 
They were predicted in numerical calculations (Ross et al. 1982) using a cylindrical model 
for magnetized plasmas with an axial current and were observed later in tokamak experi- 
ments. GAWs are excited in plasmas in a helical magnetic field, with the poloidal component 
produced by an axial current, -Bg.e/ ~ 27raBz/L, where Bz is the axial magnetic field com- 
ponent and a and L arc the plasma radius and length, respectively. The so-called tokamak 
safety factor, q = tB-^/RBq, is considered to be unity. In cold plasmas with temperatures 
less than ~ 0.1 eV, GAWs dissipate via electron- ion coUisions. However, in hot plasmas, 
when the Alfven velocity, ca — B / {Anrriiney^'^ , is of the order of the electron thermal veloc- 
ity, VTe — {Te/meY^'^, GAWs dissipatc via electron Landau damping. The electron Landau 
damping mechanism is very effective when waves interact with a group of electrons that 
have velocities approximately equal to the wave phase velocity, uj/k\\ = ca, where u is the 
frequency and fcy is the wavenumber parallel to the magnetic field lines. Electromagnetic 
field components related to GAWs are proportional to exp[i{m6 + kzZ — ujt)], with poloidal m 
and axial kz wave numbers. The requirement for these waves to propagate is that m = —1 
and kz < 0, in the case of current flowing along the magnetic field. A GAW has a very 
large quality factor, ratio of frequency to dissipation decrement, Q = cu/j ':$> 1 and a 
very long dissipation length in comparison to their wave length. La ^ Thus it is 

a very good candidate to explain plasma heating far from a source, such as in solar loops 
(HoUweg 1990; Elfimov et al. 1996). In addition a GAW is able to drive the current, which 
can be calculated by balancing the wave momentum transfer force with the electron-ion 
friction force (Elfimov et al. 1996). 
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2. Magnetosphere Accretion Model 



We study here a two-fluid magnetohydrodynamic model of a young star, accreting plasma 
(see Fig. 1). The standard parameters of CTTSs, similar to those used by Vasconcelos 
et al. (2002), are assumed: star radius i?* = 10^^ cm, truncated radius Rtr = 2 x 10^^ 
cm, truncated depth A^^ ~ 3 x 10^° cm, magnetic field at the star's surface Est = 1 kG 
and Bij. 0.1 kG in the truncated region, where the electron density is fitr — 3.0 x 10^^ 
cm~^. The length L^f of the magnetic field hues from the shock at the star's surface to the 
truncated region is ~ 1.2i?^. The velocity of the shock is Vsh ~ 2.5 x 10^ cm at the star's 
surface and Vn- ~ 3.5 x 10^ cm s~^ in the truncated region. The ionization rate parameter 

< (jv >ion for Te = 1.5 eV is ~ 3 X 10^^^ cm^ s~^ (Delcroix 1965). The accretion energy 
density is deflned by the ratio of the energy flux density to the path length L^f- Balancing 
ninety percent of the accretion power density rriinQvl^/ L^f with the ionization power density 

< av >ion ngUoEion, we obtain a plasma with an average density along the magnetic field of 
n 3.3 X 10^^ cm-^ 

We assume that the magnetic fleld connecting the truncated disk and the star's surface 
consists of filaments, along which the GAW propagates. The filament current can be driven 
by the wave ponderomotive force (see discussion below) and/or by a magnetodynamo in the 
truncated region of the disk. The current, I[a] — BaBg^ef, and current density, joh — I/'ko^ ~ 
16 statA cm^, in a filament of radius a can be produced by the electric field, j/a 10~^^ 
statV cm"-*^, that is sufficient to create the required poloidal magnetic field -Bg,e/ for the 
GAW excitation. This electric field can be generated by the induction of a perpendicular 
electric field, E± = BtrVtr/c ^ 0.05 statV cm^^, in the truncated region of the disk due 
to the magnetodynamo v x B effect (see Fig. 1). A difference of potential on the scale of 
the filament radius is of the order of 2.5 x 10^ statV. This potential is much higher than is 
necessary and can generate a filament current with a hot electron tail and a characteristic 
velocity higher than that of the Alfven velocity ca, as was found in the plasma focus device. 
We expect that the electrons in the tail will remain hot in the filament by their interaction 
with the GAW. 

For the parameters of the star assumed above, the phase velocity of a GAW, CA{Rtr), 
is ~ 3 X 10^ cm s~^ at the center of the filament. Taking into account the plasma resistance 
(see, e.g. Ginsbrug 1961), the dissipation length of the Alfven waves can be estimated to be 



La ~ Stt 



acA 



1.2 X 10^^ cm. 



(1) 



where 



(7 = 




1.6 X 10 



13^^3/2 ^-1 



= 3 X 10"^ cm-\ 
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This A; 1 1 is a maximum value, determined by the condition that the dissipation length be 
sufficient so as to explain the plasma heating along the magnetic field from the truncated 
region to the star's surface. 

Assuming a conversion of 5% of the total accreting plasma energy fiux, S — (l/2)n„mj(i'f^+ 
vf^), into the Alfven wave fiux, Sa — Pa La, we obtain a wave dissipation rate of Pa — 0.0039 

erg cm~^ s~^. This dissipation rate produces a wave momentum transfer force PAk\\/oJ (see 
Elfimov et al. 1996). Balancing the wave momentum transfer force with the electron-ion 
friction force i/gi^e^eK, we obtain the density of the driven current. 

Using a reduced friction, fei,e/ ~ ^ei{vTe/ ca)'^^'^ ■, for the electron tail with effective velocity 
~ ca and u = fcyc^ = 9 x 10^ s~^ (c^ = 3. x 10^ cm s— 1), we obtain the required current 
density, < jcd >~ 16 statA cm~^, in the filament. In this case, the electron tail that 
transports the current along the filament can be supported by the GAW via the electron 
Landau damping, while the current in the filament can be driven by the GAW. 

Here, we also present the results of calculations for the GAW in the filament. In order 
to calculate the wave field and dissipation, we use the eikonal model to analyze the wave 
propagation along the filament, exp[?'(m^^ + J kzdz — ujt)], where z is the coordinate along 
the filament axis and 6 is the poloidal angle. In this case, in order to simplify the GAW cal- 
culations using numerical codes (Amarante-Segundo et al. 1999; Galkin et al. 2002), instead 
of a dipole magnetic field, we assume a cylindrical model for the filament plasma with the 
current fiowing along a magnetic field that has a helical magnetic field line configuration. 
To calculate the wave characteristics along the filament, we take the major radius Rq — 10^^ 
cm comparable to the star's radius and the minor radius a — 5 x 10^ cm. The calculations 
are carried out with the electron density profile of the form ng — 7^0(1 ~ r'^/a^) + Ua, with 
uq = 5 X 10^^ cm~^ and = 5 x 10^° cm~^. The temperatures are taken as T^q = 1.6 
eV and Tjo = 0.5 eV. The ion density is chosen so as to satisfy the requirement of charge 
neutrality, Znz + Ui = rig, where Ui and nz are the densities of hydrogen and the other 
ions, respectively. In accordance with the Spitzer conductivity oc Te^"^ , the current density 
profile has the form j — jo(l — r'^/a'^)^. The value of the safety factor, q — rB^/RBe, is 
assumed to be unity at the filament axis and to monotonically increase up to four at the 
plasma boundary. In Fig. 2, we show the characteristic profile of the magnetic field as well 
as the dissipated power of the GAW for the filament parameters discussed above. 

In Fig. 2, we observe the peaking of the dissipation at the filament center, which predicts 
filament center peaking of the electron temperature distribution, as well. The GAW heats 
the core of the filament to higher temperature, while reducing the plasma temperature at 
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the filament border due to heat diffusion. Tfie axial wavenumber of the GAW, — uj/ca, 
varies slowly along the magnetic field lines since it depends on changes in the local magnetic 
field and density. 

3. Conclusions 

Our analysis of GAWs in the magnetospheres of CTTSs shows that: 

a) The energy of the shock (and/or accretion disk region) can be transported from the surface 
of the star (and/or accretion disk region) to the magnetosphere region by GAWs via current 

carrying filaments; 

b) The dissipation length of GAWs is on the order of the truncation radius; 

c) The conversion of only a few percent of the shock energy into GAWs is required in order 
to heat the filaments up to the observed temperature of the star's magnetosphere; 

d) The current in the filament may be induced by the electric field, created by the v x B- 
magnetodynamo effect in the truncated region of the disk. However, the major part of the 
current in the filaments is driven by GAWs. 
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VJP and RO thank the Brazihan agency GNPq and the project PRONEX (41.96.0908.00) 
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Fig. 1. — Sketch of a young accreting star of radius i?^ and truncation depth A^^- The 
electric field, designated by ©, is generated by the v x B-effect due to the interaction of the 
magnetic field lines from the star with the accreting plasma in the truncated region of the 
disk. 

Fig. 2. — Plot of Sj.^5)-components of the GAW field and dissipated power across the filament 
cross-section for a frequency of 1.4 kHz and the parameters B — IkG, qo — 1.0, no — 5x 10^^ 
cm-3, Teo = 1.6 eV, and Tjo = 0.5 eV. 
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